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Dynamical equations for time-ordered Green's 
functions: from the Keldysh time-loop contour to 
equilibrium at finite and zero temperature 

H. Ness, L. K. Dash 

Department of Physics, University of York, Heslington, York YO10 5DD, UK 
European Theoretical Spectroscopy Facility (ETSF) 

Abstract. We study the dynamical equation of the time-ordered Green's function 
at finite temperature. We show that the time-ordered Green's function obeys a 
conventional Dyson equation only at equilibrium and in the limit of zero-temperature. 
In all other cases, i.e. finite-temperature at equilibrium or non-equilibrium, the time- 
ordered Green's function obeys instead a modified Dyson equation. The derivation 
of this result is obtained from the general formalism of the non-equilibrium Green's 
functions on the Keldysh time-loop contour. At equilibrium, our result is fully 
consistent with the Matsubara temperature Green's function formalism and also 
justifies rigorously the correction terms introduced in an ad hoc way with Hedin 
and Lundqvist. Our results show that one should use the appropriate dynamical 
equation for the time-ordered Green's function when working beyond the equilibrium 
zero-temperature limit. 



PACS numbers: 71.10.Ca, 71.10.W, 71.15.m 
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1. Context 

Equilibrium, zero-temperature and finite-temperature Green's functions (GF) tech- 
niques based on many-body perturbation theory are widely used in electronic-structure 
and total-energy calculations. They are central for the calculation of the thermodynam- 
ical properties of many-body systems, as well as for the calculation of linear responses 
of systems under small time- dependent (or not) perturbations [Il[2i[3j|IJ|5j|6j[7J[E]- 

Conventional equilibrium many-body GF approaches, at zero temperature, 
concentrate on the time-ordered GF and determine its dynamical equation (in the form 
of an equation of motion, or its integral form known as the Dyson equation) using a self- 
energy which includes all the many-body interaction effects. The Matsubara formalism 
using time-ordered GF with imaginary time arguments (or imaginary Matsubara 
frequency in the dual Fourier representation) has been developed to deal with the finite 
temperature condition. Interestingly, not much has been done at finite-temperature for 
the time-ordered GF in real time or real frequency representation. 

However, when one wants to study systems driven out of equilibrium by an external 
"force" , as for example in the context of electron or heat current flow, or for any system 
driven by an external electromagnetic field (whether time- dependent or not), one needs 
to extend the equations for the dynamics of the quantum many-body interacting system 
to the non-equilibrium conditions. 

For this, the non-equilibrium Green's function (NEGF) technique [2| |9| [10 | [TT | [T2 | [T3 | 
[H] has been widely used to calculate non-linear charge and heat transport properties 
of solids, mesoscopic and nanoscopic systems. Also known as the closed time-path 
formalism [151 ED]) the NEGF formalism relies on an "artificial" time parameter that 
runs on a mathematically convenient time-loop contour (plus eventually an imaginary 
time for taking into account the initial correlation and statistical boundary conditions). 
It is a formal procedure that only has a direct physical meaning when one projects 
back the time parameters of the time-loop contour onto real times. It was introduced 
because it allows one to obtain self-consistent Dyson-like equations for the Keldysh 
Green's function using Schwinger's functional derivative technique. Transforming the 
Dyson equation to real time by varying the Keldysh time parameter over the time-loop 
contour results in a set of self-consistent equations for the different non-equilibrium 
Green's functions (advanced/retarded or lesser/greater). The NEGF technique is 
general and can treat non-equilibrium as well as equilibrium conditions, at zero and 
finite temperatures, within a single framework. 

In this paper we focus on the dynamical equation of the time-ordered GF with real 
time argument, and we show that the time-ordered GF obeys a conventional Dyson 
equation only in a limiting case: at equilibrium and zero temperature. In any other 
cases, i.e. equilibrium at finite temperature or non-equilibrium, the time-ordered GF 
obeys instead a modified Dyson equation containing correction terms. These correction 
terms were introduced in an ad hoc way in the seminal work of Hedin and Lundqvist 
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-oo+ S branch (+) 

-oo_ , ) branch (-) 
| 1 1 time 

h h t 

Figure 1. The time- loop Keldysh contour Ck for an "artificial" time parameter 
ordering on a forward time evolution (+) branch followed by a backward time evolution 
(— ) branch, the "turning" point t is taken as t — ¥ oo. 



We derive the modified Dyson equation for the time-ordered GF in a natural way 
from the Keldysh time-loop contour approach. We also show that for the equilibrium 
case at finite temperature, such correction terms can be obtained by a proper treatment 
of the imaginary time Matsubara finite-temperature Green's functions. 

The paper is organized as follows: in section El we derive the modified Dyson 
equation for the time-ordered GF from the time-loop contour formalism. We then 
concentrate on the equilibrium condition. Section |3]shows explicitly the relation between 
Hedin and Lundqvist's seminal work on real-time GFs at finite temperature and our 
result. In section HJ we show that our result can also be obtained in by using Matsubara 
temperature GFs and a proper treatment of the analytical continuation rules. We 
conclude and discuss the importance of our result in section [5j |Appendix A| and 



Appendix B contain important and useful relations for the GFs and self-energies on 



the time-loop contour. In Appendix C we provide a concrete example of our general 



results for a specific system: the single site Hubbard model at finite temperature. 

2. Non-equilibrium Green's functions on the Keldysh contour Ck 

The Green's function on the Keldysh contour Ck is defined as 

G(l,2) = -1(7^(1)^(2)), (1) 

where (1,2) stands for a composite index for space X12 and time t\$, on the time-loop 
contour Ck- The time-loop contour Ck consists of a chronological time evolution (+) 
branch followed by an anti-chronological time evolution (— ) branch [91 [TOj [111 as 
shown in Figure [U The time ordering Tc K of the product of fermion creation (w) and 
annihilation quantum fields is performed on the time-loop contour Ck- 

The Green's function obeys the equation of motion on the contour Ck [101 EE] : 

[id n - h (l)]G(l, 1') = 6(1 -!') + / d 3S(l, 3)G(3, 1'), (2) 

and the adjoint equation reads 

[-id T{ - h (l')]G(l, 1') = 6(1 -!') + / d2G(l, 2)S(2, 1'), (3) 

where ho(l) is the non-interacting Hamiltonian. 
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Hence, in a integral form, G(l, V) follows a Dyson-like equation 

G(l, 2) = g(l, 2) + J c!3d4 0(1, 3)E(3, 4)G(4, 2), (4) 

where g(l,2) is the Green's function for the non-interacting system, i.e. [id T1 — 
ho(l)]g(l, 2) = 5(1 — 2), and the time integrations run on Ck- Eq. (jlj) is valid 
provided that, in the remote past, the system is unperturbed and non-interacting and 
that the adiabatic assumption for switching on the interaction is fulfilled, as usually 
assumed within the Keldysh formalism [9|. Furthermore the non- interacting Green's 
function g(l,2) satisfies the Kubo-Martin-Schwinger boundary conditions (for example 
see Eq. (3.16) in PI). 

Upon the position of the time arguments (ti, t 2 ) on Ck , one obtains four different 
projections onto the real (normal) time arguments (ti,£ 2 ), i-e. the four "Keldysh 
components" of G(l, 1') : G f?1 ' 72 (xiti, x 2 t 2 ) with ti on the T)i — ± branch of Ck- 

So G ++ (t 1 ,t 2 ) = -i(T^'(l)^ t (2)) is the conventional time-ordered GF (though 
in the non-equilibrium conditions), and G~ + (ti, i 2 ) = — i(\I/(l)\I/^(2)) = (7(1,2) and 
(7~(£i,£ 2 ) = i( 1 I r ^(2) 1 i r (l)) = (7(1,2) are automatically ordered knowing that any 
time on the branch (— ) is later than any time on the branch (+) on Ck- Note that 
the equation of motion for (7^(1,2) does not contain the source term 5(1 — 2) since 
the time T\ and r 2 are on different branches, and in the non-interacting case, one has 
[id n - Vl)]^(l,2) = 0, and [-id T2 - /i (2)]^(l,2) = 0. 



Using the relationships between the different GFs given in Appendix |Appendix B 
one can calculate the integral equations for the different Green's functions [21 [9j [TQ1 EU 
IT2"| [T3"| [T4] . and one finds that the retarded and advanced G T,a Green's functions obey a 
Dyson-like equation: 

(7(1,2) = <f(l,2) + J d(34) /(l ) 3)E r (3,4)G r (4,2), (5) 

(likewise for G a ) where the time integrations are on the real times f^°° dtj. 

The lesser and greater G^ GFs obey the following quantum kinetic equation: 

G^(l, 2) = J d3d4 (1 + G r S r )^(l + £ a (7) + G r S^G a , (6) 

given in a compact notation. Eq.([6]) is valid within the conventional assumption 
of the Keldysh time-loop contour formalism where the interactions are switched on 
adiabatically. The most general solution for the lesser /greater Green's function can be 
found in Appendix A of Ref. [TTj. 

2.1. The time- ordered Green's function on Ck 

Following the rules of decomposition of the time-loop contour Ck onto real-time given in 



Appendix Appendix A one can expand the equation for the time-ordered GF (7 = (7 
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as follows (using compact notation) 



hence we get 



G++ = g ++ + (gZG) ++ 

= g ++ + (7) 
- g ++ Z < G > - g<Y>G ++ + 0<£" G ++ , 



G? =g l (1 + S*G" - S < G > ) 



After further manipulation, using Appendix A and Appendix B , we can show that the 



ion 



second line in Eq. (JS]) becomes ^S*G > — S > G*j = (? < (SG) > . Hence the full equati 
for the time-ordered GF is 

G t = g t (l + Y, t G t -X < G > )+g < (ZG) > . (9) 

From the equation of motion of GF on the Keldysh contour, we find 

g < (ZG) > = [ d2g<(l,2)[id t2 -h (2)]G > (2,3) , (10) 



and by using the (adjoint) equation of motion of the non- interacting GF g K , i.e. 
[— id t , 2 — h (2)]g > (1, 2) = 0, we can integrate Eq. ffTUj) by part and find that 

[^ < (EG) > ](1,3) = J dx 2 i{g < (l,x 2 t)G > (x 2 t,3)} t=+O0 -i{g < (l,x 2 t)G > (x 2 t,3)} t= ^ oo , 

(11) 

which represents a boundary condition term for the time variable (after all we are solving 
a first order differential equation in time). 
Hence Eq. (jHJ) is transformed into 

G*(l,2) = < 7 t (l,2) + f d(34) g\l,3) (S'(3,4)G*(4,2)-S<(3,4)G > (4,2)) 

J (12) 
- / dx 2 i[g < (l,x 2 t)G > (x 2 t,3)] t=+oo -i[g < (l,x 2 t)G > (x 2 t,3)] t= _ O0 , 

or equivalently (in compact notation) 

G f = g* (l + - S<G>) + boundary condition. (13) 

Eq. (ED, equivalently Eq. OH]), is the most general result for the time-ordered Green's 
function and is valid for the non-equilibrium case at finite and zero temperature. It also 
clearly shows that the time-ordered Green's function G*(l,2) does not follow a Dyson- 
like equation. However this is not a closed equation, in contradistinction with the 
Dyson equations for the advanced and retarded components, since it also involves greater 
and lesser components for the GF and the self-energy. In the general non-equilibrium 
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conditions, one rather works with the closed equations obtained for the retarded and 
advanced GF which are related to spectroscopic properties, and with the quantum 
kinetic equations obtained for the lesser and greater GF. Though at equilibrium, the 
lesser and greater components of GF or a self-energy are related to the respective spectral 
functions and the equilibrium distribution. Hence we recover closed equations for the 
time-ordered GF. 

We also establish in Section H] the formal links with the approach of Hedin and 
Lundqvist and with the Matsubara temperature GF formalism. 

It would be interesting to find some conditions under which the the boundary 
condition term in Eq. ( fT3l) vanishes. At equilibrium, because of the principle of adiabatic 
switching (on/off) of the interactions in the remote past/future, the initial state returns 
(up to a phase factor) to itself in the far future. Therefore it would be tempting to 
consider that the the two terms defining the boundary condition are the same at t = ±00 
and that the boundary condition term vanishes. However this is not that shown 
in Section HI 

At non-equilibrium, the problem is different and the solution may depend on the 
size of the system. For an infinite system (a large system in the thermodynamic limit), 
in some cases, we would like to use the so-called memory- loss theorem [T7J [181 US]: the 
long time limit of the lesser (or greater) components of the GF or of the self-energy 
vanish. Hence in that case, the boundary condition term may vanish (though this 
needs to be proven). This corresponds to cases when the system "thermalise with the 
environment" . For finite size systems at non-equilibrium, no "thermalisation" is possible 
and the boundary condition term play an important role in the solution of dynamical 
equation of the time-ordered GF. 

Before turning onto the equilibrium condition, let us make two further general 
remarks. One could note that from the definition G r (l,2) = (7* (1,2) — G < (1,2) 
(see Appendix Appendix B ) or equivalently G t (l,2) = G r (l,2) + G < (1,2), one could 
have gueseds that G l does not obey a Dyson equation since G r does, but not G < 
|20j . Futhermore we have also checked that we can also recover, as expected, our 
result Eq. ffl3|) by using the Kadanoff-Baym formalism which is somewhat different but 
equivalent to the Keldysh time-loop approach [2Tj . 



2.2. At equilibrium 

We now turn onto the special case of equilibrium for which the Green's functions 
techniques based on many-body perturbation theory were orginally derived. 

At equilibriumn, as in the steady state, all quantities X(t, t') (GFs and self-energies) 
depend only on the time difference X(t—t'). One can then Fourier transform the quantity 
X(t — t') and work in the real- frequency representation X(u). 

The main result of the paper, Eq. ( Tl2|) . is recast as: 

G\u) = g\u) (1 + S*(w)G*(a;) - - 0<(w)(E(w)G(w))> , (14) 
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in compact notation (i.e. omitting the explicit integrals over space). 

It should be noticed that the last term S < G > in Eq. (j!2p satisfies the detailed 
balance equation at equilibrium [TO]: H < (oj)G > (u)) = T l > (u)G < (u). 

Eq. (fT2l) can also be recast as 

[[^(c,)]- 1 - E»] G*(w) + E<( W )G>( W ) = 1 , (15) 

ignoring the boundary conditions term. 

Furthermore, at equilibrium, the lesser and greater components of either a GF or a 
self-energy (X^) can be expressed in terms of the corresponding advanced and retarded 
quantity and a distribution function [SI SI E21 El [231 El E31 ESI ES] , i.e. 

X*{u) = -f^(u)(X r (u) - X a (u)), (16) 

and spectral function A x = {X r {u) - X a (w))/(-27ri). 

At equilibrium and for a system of fermions, / <,eq is given by the Fermi-Dirac 
distribution function f c %u) = [1+exp ^(uj-^)}^ 1 and />> eq = / cq -l (with (3 = 1/kT). 
Eq. ( flBT) follows from the Kubo-Martin-Schwinger relationship between the lesser and 
greater component at equilibrium: X >,cq (u;) = — e f3 ^^^X < ' cq (uj). 

Hence we can easily see that 

E<(u;)G» = (27r) 2 r q H(l - r(u))A s (u;)A G (u;). (17) 

Therefore in the limit of zero temperature, there is no contribution from £ < G > 
since / cq (l — / cq ) = (this property holds for any product X < Y > , and hence for the 
boundary condition term g < (EG) > ), and the time-ordered GF G l obeys a conventional 
Dyson-like equation G\u) = g\uj)[l + Y>\uj)G\uj)]. 

It is only in the case of equilibrium at zero temperature that G l follows a Dyson 
equation otherwise at finite temperature G t obeys a modified Dyson equation, Eq. (ITS]) , 
including the extra term E < G > and the boundary condition term g < (EG) > . 

In the following section, we are going to show how Hedin and Lundqvist noticed 
the necessity to introduce correction terms in the real-time Green's function at finite 
temperature [16]. They found it necessary to introduce extra terms in the Dyson 
equation of the real-time Green's function at finite temperature to recover the expected 
behaviour of the independent particle case. The terms introduced ad hoc by Hedin and 
Lundqvist correspond exactly to the E < G > term that we have obtained in a natural 
way from the use of the time-loop contour. 

3. Real-time Green functions at finite temperature 

In this section, we retrace the analysis made by Hedin and Lundqvist (section 17 of 
Ref. [16]) on the temperature Green's functions. The definition for the real-time GF at 
finite temperature follows from the zero temperature case by using an ensemble average 
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for the time-ordered product of fermion operators instead of a ground state expectation 
value, i.e. 

G(l,2) = -iTr[p G T*(l)¥(2)] , (18) 

where the statistical operator (density matrix) p G = Zq 1 exp — (3(H — p cq N) and the 
partition function Z G = Tr[pc] are given for the grand canonical ensemble. 

After Fourier transform, we can write the time-ordered GF using a spectral 
representation just as in the zero temperature case [3j HJ 116] : 



G{oj) = / doj'A(oj') 



(19) 



oj — oj 1 + ir] oj — oj 1 — ir] 

where we omit the spatial variables. 

From Eq. (fT9l) . we can obtain the real and imaginary parts of G{oj): 

^eG{oj) = PV [ doj'^^, (20) 
J oj-oj' 

and 

3mG(u) = n[2f cq (u) - l]A(u) . (21) 

Hedin and Lundqvist noticed that by considering the conventional Dyson equation 
for the time-ordered GF: [oj — e — H(oj)]G(oj) = 1, one derives the following spectral 
function 

B(oj) 

^ ~ [oj-e- Ke£(u;)] 2 + Ti 2 B(ojf[2f^ - l] 2 ' ^ 

with the spectral representation for the self-energy SmE(w) = 7i[2f cq (oj) — l]B(u). 

However Eq. ( 1221) does not have the proper limit in 5(oj — e) for the independent 
particle case where 3?eS(o;) = and B(oj) — > 0. 

To palliate this discrepancy, Hedin and Lundqvist introduced ad hoc a correction 
term in the conventional Dyson equation that G and S must obey. Their modified 
Dyson equation reads 

[oj-e- £(w)]G(w) + (27r) 2 / cq (l - f C(l )A(oj)B(oj) = 1. (23) 

From this modified Dyson equation, one can calculate the imaginary part of G{oj) and 
obtain the spectral function A(oj) 

A ^ = r » ■ 2u ( yi > (24) 

[oj — e — UieL(oj)\ z + tx 1 B{oj) 1 

as for the zero temperature limit, it has the proper limit Aioj) — > 5(oj — e) for the limit 
KeS(w) = and B(oj) -> 0. 

Hedin and Lundqvist conclude "that the rather messy definition of S in Eq. ( d3[) 
shows that it is better if, at finite temperature, we can work with concepts other that 
the straight-forward generalisation in Eq. (TT2j] of the zero temperature Green's function. 
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One possibility is to use the Matsubara technique with temperature Green's functions, or 
the close related Martin and Schwinger technique ..." [T6] . 

It is interesting to note that the "messy" definition Eq. ( 123]) is actually the correct 
result for the time-ordered GF G l at finite temperature. Eq. (I23j) is strictly equivalent 
to Eq. ( Fl5l) when one notices that in the definition of S < G > in Eq. ( fT71) one deals with 
the same spectral functions, i.e. A = A G and B = A 2 . 

So by using the Keldysh time-loop contour, we have found in a natural way, almost 
straightforward, the proper modified Dyson equation for G l at equilibrium and this, 
without the need of introducing correction terms with no proper justification, other 
than the need to recover the correct behaviour for A G (u) in the non- interacting case. 

To conclude this section, one could rephrase our result by saying that the time- 
ordered GF G l at equilibrium and finite temperature obeys a Dyson equation G l = 
g l + g t 'EG t only if one defines the corresponding self-energy £ as E = E* — H < G > (G t )~ 1 . 

4. Recovering G t from the equilibrium formalism and imaginary-time 
temperature Green's functions 

In the next section, we show that the proper modified Dyson equation for G f can be 
obtained from the knowledge of the spectral function A x {uj) = (X r (u) — X a (u))/ — i2n, 
with X being a GF or a self-energy. We also discuss the connection of our result with 
the Matsubara temperature Green's functions. 

We extract the spectral function A G using the Dyson equations for G r ' a : 



Using the definitions of the spectral functions A G and A 2 and the fact that 
^eG r = $leG a = P.V. J A g (lo')/(lo - u') = fteG* as in Eq. ([20]) (with similar relations 
for the self-energy £), we find 



The time-ordered GF G f (u) can then be obtained from its spectral representation 
Eq. ffT9]) using the expression of the spectral function A G {u) given by Eq. f )26|) . The 
first term in the RHS of Eq. ( 126]) generates the non-interacting time-ordered GF 
g*(o;), while the two other terms regrouped in D(u) = A 9 ($leY>MeG — ^A^A ) + 
3fte(g).(9?eE.A G + 9fteG.A s ) give a real part P.V. f D(uj')/(uj — uj 1 ) and an imaginary 
part 7r[2/ cq (u;) — l]D(u). The proper result for the modified Dyson equation of G*(u;) 
can be obtained after lengthy algebraic derivations. Instead we now show another just 
as valid shorter route to end up with the correct result. 

We consider the different contributions of a conventional Dyson g l + g t H t G t . From 
the spectral representation, we have already generated the term in g f from A 9 , so we 
are left with the g^G 1 . 



G r -G a = g r -g a + g r YI'G r - g a E a G a . 



(25) 



A 9 + A 9 (jR,eT,MeG - n 2 A*A G ) 
+ ^e(g). (jReY,.A G + A^MeG) . 



(26) 
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= (KeS + i7r(2/ cq - 1)A S ) (jfteG + i7r(2/ eq - l)A G ) 



(27) 



= (KeS.KeG - n 2 A^A G - (27r) 2 / eq (/ eq - l)A^A G ) 

+ ivr(2/ eq - 1) (9tfe£.A G + A^MeG) 
= (KeS.KeG - ti 2 A^A g + S < G > ) 

+ ivr(2/ eq - 1) (9tfe£.A G + A E 9fJeG) , 

where we have introduced, in the last equality, the definition of the product T, < G > given 
in Eq. (I17j) . We can see that the correction term in £ < G > appears already at this level. 
Hence one can see that 

1 + S*G* - S < G > 

f 28) 

= (1 + IkeTMeG - vr 2 A s A G ) + i(2/ oq - l)vr (jR,eY,.A G + A^eG) . 

To finish our derivation, we first consider the imaginary part of g*(l + £*G* — S < G > ). 
By using the definition for g f from Eqs. ( )20l) and ( 12TT) = 9?e(g) + i(2f eq — l)irA 9 } we 
obtain: 



3m [g\l + E*G* - E<G ; 
= 3m(0*) 9fce(l + E'G* - S < G > ) + Ke(^) Qm(l + E'G" - S < G > ) 
= (2/ oq - l)nA 9 (1 + MemeG - tc 2 A e A g ) + Ke(#)(2/ cq - 1)tt (3?eS.A G + A s 9fceG) 
= (2/ cq - l)ixA G 
= SmG 1 . 

(29) 

For the real part of g\l + E t G t - £<£>), we find that 



9f?e + E'G* - E<G° 
= Ke(#*) + E*G* - S < G > ) - 3m( 5 ( ) Qm(l + E'G* - S < G > ) 
= Ke(^) (1 + VlemeG - n 2 A*A G ) - (2/ cq - l)vrA 3 (2/ cq - 1)tt (KeE.A G + A^eG) 
= Me(g) (1 + IkeTMeG - tt 2 A s A g ) - tt 2 A 9 (9£eE.A G + A s KeG) 
_ 4/ cq (/ cq - 1)txAH (3?eE.A G + A^eG) . 

(30) 

We can use the properties of the Hilbert transform [27JEEJI29] , i.e. %[H[f] ] = —f and 
the Bedrosian indentity %[fg\\ = fH[g] = gH[f], to show that: 

ReG* = U[ixA G } 

= H\tiA 9 (1 + IfreYMeG - tt 2 A^A g )} + H[irUe(g) (jRe'E.A G + A s 3?eG)] 

= H[nA 9 } (1 + KeEKeG - n 2 A s A G ) + U[n$le(g)} (9fJeE.A G + A^eG)} ' ° 

= tte( ff ) (1 + 9£eE9£eG - vr 2 A E A G ) - vr 2 A 9 (9£eE.A G + A^eG)} , 
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since &e(g) = n[nA9], and H[n$le(g)] = k 2 U[U[A9] ]. 
Hence Eq. ( 151]) becomes 

9fJeG" = 3fJe [#*(1 + S*G" - S < G > )] + 4/ cq (/ cq - l)nAH (KeS.A G + A s KeG) , (32) 

When considering that, by definition at equilibrium, we have: 

g < = _ r{g r _ ga) = ir27rA9 

(EG) > = -(/ cq - 1) ((SG) r - (SG) a ) , 

and (£G) r - (SG) a = S r G r - S a G a = -2i{n A^eG + Ke£.7rA G ), we can see that, as 
expect, the last term in Eq. (|32|) is just equal to g < (EG) > . 

Hence we recover once more the modified Dyson equation for the time-ordered GF: 

G\u) = g\u) (1 + E'MG'M " S < (w)G > (a;)) - 9 < (uj){Y,(uj)G{uj)) > (34) 

in a indirect way from the spectral functions. 

We now briefly comment on the imaginary-time temperature Green's functions. 
The definition of the temperature GF is, in form, similar to a time-ordered Green's 
function: 

£(xr, xV) = -Tr [p G T r ^(xr)^t(xV)] , (35) 

where one uses imaginary times t = it instead of real time t, and the T T orders the 
operators according to their value of r, with the smallest at the right [H H] . 

The temperature Green's function Q(l,2) follows a conventional Dyson-like 
equation: 

0(1, 2) = O (1, 2) + J d(34) g (l, 3)S(3, 4)£(4, 2), (36) 

where the time integrations run over Tj from to /3 (and Qq is the temperature Green's 
function for the non-interacting system). 

At equilibrium, the quantities depend only on the "time" difference T\ — r 2 and 
(?(x, x',r — t') is periodic over the range 2/3 (or periodic (antiperiodic) over the range 
[0,/3] for boson (fermion)). Therefore Q can be expanded in a Fourier series with 
coefficients C/(x, x', cj n ) given in terms of the Matsubara frequencies u n being integer 
even (odd) numbers of 7r//3 for boson (fermion). Each Fourier component obeys also a 
Dyson equation: 

G(u n ) = £o(w n ) + £ (w„)£(w n )£(w n ) (37) 

(in compact notation). 

There is a well defined mathematical connection between the temperature GF Q(u n ) 
and the real-time GF, this is obtained from the analytical continuation from the complex 
frequency iu n on the real frequency u (and changing, when needed, the sum over the 
Matsubara frequencies onto integrals 1/(3 •••—)■ J dcu/27i . . .). 

However, even if Q is an imaginary-time-ordered GF, it does not transform into 
the real-time-ordered GF. The analytical continuation iu n — » u + irj transforms the 
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temperature GF into the retarded real-time GF (in energy representation) [21 H] : 
Q{oj n ) G r {to). Similarly we have Q{u n ) — > G a {oj) for iuj n — > u — ir]. This is also 
consistent with the fact the G r,a obeys a Dyson equation, as much as Q does too. 

In order to get the time-ordered GF G*(a;) for the temperature GF Q(u n ), one has 
to adopt the following procedure [I]: from Q{u n ) and the analytical continuation, one 
obtains G t (uj) and G a (u) from which one extracts the spectral function A (u) which is 
then used to built the time-ordered GF from the spectral representation Eq. ffl9|) . 

Following this procedure, we have shown above that one can obtain the modified Dyson 
equation for G"(w). 

As another example, we provide in Appendix C similar calculations for a specific 
system, i.e. a finite size system at finite temperature, described by the Hubbard 
Hamiltonian. 



5. Conclusion 

We have derived the dynamical equation for the time-ordered GF in the most 

general case of non-equilibrium and finite temperature. In the equilibrium case, we have 
shown that the corresponding time-ordered GF G t (u) obeys a Dyson equation only in 
the limit of zero temperature. At finite temperature, the Dyson equation for G 1 {uj) needs 
to be modified by a correction term. We have shown that the correction term correspond 
exactly to the corrections introduced ad hoc by Hedin and Lundqvist in their treatment 
of the real-time GF at finite temperature. We have also shown that these corrections can 
be obtained from the Matsubara temperature GF following the appropriate protocol for 
the analytical continuation from the imaginary frequency onto the real frequency and 
the spectral representation of the time-ordered GF. Table [T] summarizes the different 
GFs formalism and the corresponding dynamical equations for the different GFs. 

Our work leads to two important conclusions. First, the Keldysh time- loop contour 
formalism to the GFs, though developed to deal with non-equilibrium conditions, is 
also extremely appropriate to the study of the equilibrium properties of the system. It 
leads, in a more natural and straightforward way, to the correct dynamical equations 
of the GFs. Second, any formalisms based on the resolution of the Dyson equation for 
the time-ordered GF, as those used for example in electronic structure and excitations 
calculations, are valid only at equilibrium and in the zero temperature limit. These 
formalisms include GF calculations ranging from the GW approximation [161 EH ED] to 
more recent work on self-consistent GF equations involving four-point propagators [32J . 
Beyond the zero temperature, correction terms in the Dyson equation for G l (u) must 
be taken into account, or one needs to work with the Matsubara formalism for Q(iu n ). 
Such temperature effects were also found for the T-matrix approximation in a real-time 
formalism in Ref. 1331. 
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Non-equilibrium 
all T 



T ^ 

at equilibrium 



T = 

at equilibrium 



G(t,t>) 

G r ,G a 
in (t, t') or 

in (£, £') or 

( W ) 

67* 

in (£,£') or 



in (ir, ir') or 

(it*; n ) 



Dyson eqn on Ck '■ 
G = g + gZG 

Dyson eqn: 

Qr,a _ gr,a _|_ gr,aj^r,aQr,a 



QKE: 

+(l + G r £ r )^(l + £ a G' a ) 

Modified Dyson eqn: 

G l = ^(1+S*G*-S < G > ) 

+ boundary condition 



Dyson eqn: 

Qr,a _ gT,a _|_ gr,aj^r,aQr,a 



QKE: 

G^(w) = 

-f^(u)(G r - G a )(u) 

Modified Dyson eqn: 

G l = ^(l+S*G t -S < G > ) 

+g<(ZG)> 

Dyson eqn: 

Q = Go + Go^Q 



Dyson eqn: 

Qr,a = 

gT,a _|_ gr,aj^r,aQr,a 

QKE: 

G* = 

Dyson eqn: 
G l = g l + 



Table 1. Summary of the different Green's function formalisms in different situations 
(equilibrium/non-equilibrium, finite/zero temperature). The rows show the relevant 
Green's functions and the different dynamical equations (Dyson, modified Dyson, 
or Quantum kinetic equation QKE) for each situation. Compact notation is used, 
involving integration of space and time (unless the energy uj representation is used) 
coordinates. 
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Appendix A. Rules for continuation on time-loop contour Ck 

For the following products P(i)(t, t') on the time-loop contour Ck, 

P (2) (r,r')= / df A(t,t)B(T,t') 
Jc K 

P m = [ ABC 
Jc K 



P(n) = / A 1 A 2 ...A n , 
Jc K 



we have the following rules for the different components P%Jt,t') on the real-time axis: 

(x = r, a, >, <) 

P ( |(t, t') = J di A r (t, i)B^{t, t') + A*{t, i)B a (i, t') 
p ( | = / A^B a C a + A r B^C a + A r B r C^ 



pr / at at at pa _ ACl ACl ACl 



(n) ~ / rx \ rx 2---- rL n 1 1 (n) 

Appendix B. Relationship between the different Green's functions and 
self-energies 

The relations between the different components of the Green's functions and self- 
energies on the Keldysh time-loop contour are given as usual, with X mm (12) = 
GWi 1 »(12),E" 11 »(12) or g™{Yl). 

X r = X 1 - x< = x> - x i 

X a = X 1 - X> = X< - X 1 

B.l) 

X 1 + X 1 = x< + x> 
x> -x< = x r -x a 

The usual lesser and greater projections are defined respectively as X < = X + ~ and 
X > = X~ + , and the usual time-ordered (anti-time-ordered) as X 1 = X ++ (X f = X ). 

Appendix C. Example for a finite system at equilibrium 

In this section, we consider in finite system, with interaction, for which analytical 
expressions are derivable for the different GFs and we check the validity of Eq. (|T4"I) 
for such a system at equilibrium. 

We work with the single-site Hubbard Hamiltonian given by: 

H = ^2 £ ^ n ° + Ur H n i (C- 1 ) 

<T=t4 
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where n a is the occupation number operator of the spin cr. In the absence of magnetic 
field, the system is spin degenerate when e Qa = e Q . We consider this case in the following. 
The system has 4 different states upon the level occupation n a = 0, = 1 and nj, = 0, 
= and n± = 1, = 1 and n± = 1 with respective energy 0, e , £o, 2eo + U. One can 
calculate the advanced and retarded GFs at equilibrium (and finite temperature, for the 
finite size system "connected" to a thermal bath at temperature T) following the usual 
prescription [3JH]. Within the Lehmann representation, we find 

1 f 1 + p -/3eo p -/3eo i p -/3(2e +U) } 

G r '» = j + , . + + jjTT. , (C.2) 

z [ a; — £o ± u — e — U ± if] ) 

with the partition function Z = 1 + 2e _ ^ eo + e - /3 ( 2eo+f/ ), /3 = 1/fcT and 77 ->■ + . 
We define the thermal occupancy / as follows 

/ = - ^ and l-/ = i±|_. (C.3) 

In the non- interacting case (U = 0), one recovers the usual GF expressions for a 
single electronic level, G r,a = 1/{oj — £0 ± irj) and / = / eq = 1/(1 + e^ eo ) the Fermi 
distribution (with the chemical potential /i eq = taken as the reference of energy for all 
the expression given in this Appendix). 

Eq. (1C.2I) can be rewritten as 

G r ' a (u) = l ~ f + 1—— , (C.4) 

OJ — Eq ± 17] OJ — S — U ± 17] 

with the usual relationships between the imaginary and real of G r ' a (oj) based on the 
properties of the Hilbert transform %[k8{oj — E)\ = VV -^z^s (with W standing for the 
Cauchy principal value). 

We now define the self-energy Y7' a (oj) corresponding to G r ' a (oj) as follows: 

+ f - = (C5) 



OJ — Eo ± 17] OJ — Eo — U ± 17] OJ — Eq — H r ' a {oj) 

One can separate the real and imaginary parts of both sides of the equality, using 
E r / a = KeS ± i^mS: 

(1 - /)PV— L_ + fVV - =F ivr [(1 - f)S(uj - e ) + fS(uj - E - U)} 

U — Eo OJ — E — U , . 

oj — Eq — 3fteS =p i^mE ^ ' ' 
= {oj-Eq- KeS) 2 + (3mS) 2 ' 

In the LHS of Eq. (1C.6I) . the real and imaginary parts are related to each other via 
an Hilbert transform as explained above. Hence this must held as well for the RHS of 
Eq. flC.6j) . Indeed, from Eq. (IC.6j) . we have 

oj - sq - 9£e£ = (1 - f)VV—^— + fVV x N{oj) , (C.7) 

OJ — Eq OJ — Eo — U 
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and 

SmS = [(1 - f)5(u - s ) + fS(u - £ - U)] x Af(u) , (C.8) 
with Af-^u) = (u - e - KeS) 2 + (SmS) 2 . Combining Eqs.flGl]) and (lCl8]) r we find 



(1 - f)VV—^— + /PV- 1 



+ [(1 - - £ ) + /^(w - £o — A/" 2 

(C.9) 

hence another definition for f\f(u). 

By using three properties of the Hilbert transform [27[ EEJ [29] , T-L[H[g(u)] ] = —g(u), 
and the Bedrosian identity for the Hilbert transform of product functions H[g(u)h(u)] = 
g(u)T-L[h(u)] = hiv^T-ilgiu)}, one can show that from Eqs. (IC.7j) and (IC.8j) . we have: 

H[(u - e - 3?eS)A/] = -U[UeY] = ( (1 - f)U[VV—^— ] + /ft[PV -] J M 

V ^ — £ CO — 6 — U J 

= -SmS , 

(CIO) 

hence 9fte£(u;) = -K[3mS(w)] as expected, see for example Eq. fl20l . The first equality 
comes from the fact that %[(cj — eo)N] = when A/" _1 (w = £q) ^ 0. 

Now that we have established the usual relationships between the imaginary and 
real parts of both the advanced and retarded GFs and self-energies (X) for our finite 
size model, i.e. X r ' a (u) = ^eX(u) q= i3mX(w) and X r (u) - X a (co) = -2niA x (with 
A x (oo) = ±iSmX r / n (w)), and ^eX(uj) = VV Jdu'^ff- = -U[%mX(uj)}, we can go 
further along and define as usual the lesser and greater components (at equilibrium): 
X<(u) = -f{X r {u)-X a {u)) = 2irifA x {u),&ndX > {u) = -(/ - l){X r {u) -X a {u)) = 
2ni(f — l)A x (u), as well as the time-ordered components: 9fteX'(u;) = 9fteX r,a (o;) and 
^smX t (uj) = 7r(2/ — \)A x {uj). We recall that for our finite size system at equilibrium, 
the distribution / is given by Eq. flC.3j) . 

We can now check the validity of Eq. (|14p for our model system. For that we 
follow the prescriptions given in Section [H We first evaluate the imaginary part of 
g* + g^G 1 - g t Y < G > : 

3m [g* + g'T^G 1 - ^S < G > ] = 

3m 5 ' (1 + KeE^eG* - SmE'SmG* - S<G > ) + Ueg 1 (SmS'KeG* + KeS'SmG") , 

(C.ll) 

since H < G > is real. Using the different definitions, we find 
3m [g* + #*E*G* - g*E < G > ] = 

$W (1 + 'ReYMeG - n 2 (2f - l) 2 A^A G - (2n) 2 f(l - f)A^A G ) (C.12) 
+ Meg 1 (tt(2/ - l)A^eG + KeE7r(2/ - l)A G ) . 
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Using the definition Qmg 1 = tt(2/ — 1)A 9 , we end up with 
3m [g* + y*E*G* - ^S<G>] = 

tt(2/ - 1) [A 9 (1 + MemeG - n 2 A^A G ) + Meg (A^eG + KeSA G )] , 

where the term in brakets is just the definition of the spectral function A G given by 
Eq. ( |26l) . Hence we find that 

3m [</* + g^G 1 - ^S < G > ] = tt(2/ - 1)A G = 3mG" , (C.14) 

as expected. One can proceed along the same lines to establish the relation for the 
real parts, using the different properties of the Hilbert transform mentioned above. 
Therefore, as expected, Eq. ( !T4|) is valid for the finite size system, at equilibrium, 
described by Hamiltonian Eq. (IC.lj) . 
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